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Unresolved emission and ionized gas in the bulge of M31 
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ABSTRACT 

We study the origin of unresolved X-ray emission from the bulge of M3 1 based on archival 
Chandra and XMM-Newton observations. We demonstrate that three different components 
are present: (i) Broad-band emission from a large number of faint sources - mainly accreting 
white dwarfs and active binaries, associated with the old stellar population, similar to the 
Galactic Ridge X-ray emission of the Milky Way. The X-ray to K-band luminosity ratios 
are compatible with those for the Milky Way and for M32, in the 2-10 keV band it is 
(3.6 ± 0.2) x 10 27 erg s" 1 Lq. (ii) Soft emission from ionized gas with temperature of about 

~ 300 eV and mass of ~ 2 x 10 6 M©. The gas distribution is significantly extended along the 
minor axis of the galaxy suggesting that it may be outflowing in the direction perpendicular 
to the galactic disk. The mass and energy supply from evolved stars and type la supernovae 
is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral 
arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the "vertical" 
direction, (iii) Hard extended emission from spiral arms, most likely associated with young 
stellar objects and young stars located in the star-forming regions. The Lx/SFR ratio equals 
~ 9 x 10 38 (erg/s)/(M©/yr) which is about ~ 1/3 of the HMXBs contribution, determined 
earlier from Chandra observations of other nearby galaxies. 
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1 INTRODUCTION 

The X-ray radiation from the majority of galaxies is domi- 
nated by X-ray binaries (e.g. Fabbiano 2006). In addition, extended 
emission is present in galaxies of all morphological types. At least 
part of this emission is associated with stellar population and is a 
superposition of a large number of faint compact sources - accret- 
ing white dwarfs, active binaries and other types of stellar sources, 
old and young (Revnivtsev et al. 2006; Sazonov et al. 2006). There 
is also a truly diffuse component - emission from ionized gas of 
sub-keV temperature. Its importance varies from galaxy to galaxy, 
with luminous gas-rich ellipticals, like NGC 1316 (Kim & Fab- 
biano 2003) and dwarf galaxies similar to M32 (Revnivtsev et al. 
2007) representing the two opposite ends of the range. Signifi- 
cant diffuse or quasi-diffuse emission is also associated with star- 
formation, the Antennae (e.g. Baldi et al. 2006) being one of the 
nearby examples. The morphology of the gas in starburst galaxies 
often indicates outflows, driven by the energy input into ISM from 
core collapse supernovae. Theoretical considerations suggest that 
gas in low mass elliptical galaxies may also be in the state of out- 
flow (David et al. 2006). The mass and energy budget of the ISM 
in this case is maintained by winds from evolved stars and Type 
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la supernovae. The overall X-ray radiation from a galaxy is a su- 
perposition of these (and possibly other) components, their relative 
importance being defined by the morphological type of the galaxy 
and its star-formation and merger history. 

Our close-by neighbor, M31 galaxy gives a unique opportu- 
nity to explore a "full-size" spiral galaxy similar to the Milky Way 
without complications brought about by projection and absorption 
effects, often hampering studies of our own Galaxy. Not surpris- 
ingly, it has been extensively investigated by every major observa- 
tory of the past decades. Observations by the Einstein observatory 
demonstrated that X-ray binaries account for the most of the X-ray 
emission from the galaxy (van Speybroeck et al. 1979; Fabbiano 
et al. 1987). Using the complete set of Einstein data Trinchieri & 
Fabbiano (1991) constrained possible amount of ionized gas in the 
bulge of the galaxy by < 2 x 10 6 M o . Based on the ROSAT ob- 
servations, Primini et al. (1993) found evidence for extended emis- 
sion component with luminosity of ~ 6 X 10 38 erg/s in the 0.2- 
4 keV band. They suggested, that this emission may be of truly 
diffuse origin or due to a new class of X-ray sources. The unre- 
solved emission from M31 was investigated further by Supper et 
al. (1997), West et al. (1997), Irwin & Bregman (1999), Borozdin 
& Priedhorsky (2000). In all these studies the existence of a soft 
emission component with temperature kT ~ 0.3 - 0.4 keV has been 
confirmed, although different authors suggested different explana- 
tions of its origin. With advent of Chandra and XMM-Newton, the 
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consensus seemed to be achieved in favor of truly diffuse origin of 
the soft emission component (Shirey et al. 2001; Takahashi et al. 
2004). However, recent progress in understanding the nature of the 
Galactic Ridge emission as a superposition of a large number of 
faint stellar type sources (Revnivtsev et al. 2006), made it worth 
to revisit the problem of the origin of extended emission in M31. 
Also, with more Chandra and XMM-Newton observations, more 
accurate and detailed investigations became possible. Recently (Li 
& Wang 2007) analyzed large Chandra dataset of M3 1 bulge obser- 
vations and demonstrated presence of both ionized gas and emis- 
sion of faint compact sources associated with old stellar population. 
Moreover, they revealed peculiar morphology of the gas emission 
and suggested that the X-ray gas in the bulge of M3 1 may be in the 
state of outflow. 

In the present paper we combine extensive set of Chandra and 
XMM-Newton observations to obtain a broad band and large field 
view of the X-ray emission originating in and around the bulge of 
M31. We restrict our study to the central region of ~ 20' in radius, 
covered well by Chandra and XMM-Newton observations currently 
available in the public archives of these observatories (Fig.l). The 
investigated region has the linear size of ~ 4 kpc along the major 
axis of the galaxy, but extends out to ~ 16 kpc along the plane of 
the galaxy in the minor axis direction, due to rather large inclination 
angle of M31, i ~ 77° (Henderson 1979). We assume the distance 
to Andromeda of 780 kpc (Stanek & Garnavich 1998; Macri 1983). 
The Galactic absorption towards M31 is 6.7 x 10 20 cm -2 (Dickey 
& Lockman 1990). 

The paper is structured as follows. In Section 2 we describe 
the data and its reduction. We introduce our results in Section 3, 
where we present the spatial distribution, morphology and the spec- 
tra of the extended emission. In Section 4 the origin, properties and 
physical parameters of different components are discussed and in 
Section 5 we summarize our results. 



2 DATA REDUCTION 

We combine data from Chandra and XMM-Newton satellites 
adding their benefits together. The primal advantage of Chandra is 
its good angular resolution which allows us to resolve individual 
X-ray binaries everywhere including the very central region of the 
bulge. XMM-Newton provided better coverage of M31 and col- 
lected more photons, thanks to its larger effective area. It is more 
suitable to study the outer part of M31. On the other hand the 
higher and less predictable background of XMM-Newton compli- 
cates study of low surface brightness regions. 



2.1 Chandra 

We used 21 archival Chandra observations listed in Table 1, 
taken between 13.10.1999 and 23.05.2004. For the analysis we ex- 
tracted data of the ACIS-I array except for OBS-ID 1575 where we 
used only the S3 chip. The pattern of available Chandra observa- 
tions allows us to study the central « 15' region. The data reduc- 
tion was performed using standard CIAO 1 software package tools 
(CIAO version 3.4; CALDB version 3.4.1). For each observations 
we filtered out the flare contaminated intervals, excluding the time 
intervals where the count rate deviated by more than 20% from the 
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Figure 1. RGB image of M31. The colors are as follows: red is the XMM- 
Newton data in the 0.5 -1.2 keV band, green is an image of the K-band 
light from 2MASS and blue is the 160 /urn image of Spitzer. Center of the 
M31 is marked with a cross. North is up and east is left. 



mean value. The resulting effective exposure times are given in Ta- 
ble 1. 

Crucial for the analysis of the low surface brightness outer 
regions is the accuracy of the background subtraction. In treat- 
ing the Chandra background we generally followed the procedures 
outlined in Hickox & Markevitch (2006). We determined the level 
of the instrumental background using the stowed data set 2 . In the 
stowed position of ACIS detectors the sky emission is blocked 
and only the instrumental background gives a contribution. As 
demonstrated in the above-mentioned paper, although the instru- 
mental background level varies with time, its spectrum remains un- 
changed. The effective area of Chandra is negligible above 9 keV 
and the count rate is dominated by the instrumental background. 
Therefore the 9.5-12 keV count rates can be used to renormalize 
background spectra obtained from the stowed datasets. In dealing 
with cosmic background we took into account that it consists of the 
soft emission associated with the Galaxy and harder extragalactic 
component and treated them separately. For the soft Galactic com- 
ponent we used the best fit spectrum from Hickox & Markevitch 
(2006). For the extragalactic background we took into account that 
some fraction of it has been resolved in our data analysis procedure 
and removed along with X-ray binaries. Using the incompleteness 
function from Voss & Gilfanov (2007), that was obtained using es- 
sentially the same data set, we estimated that our point source de- 
tection sensitivity in the outer regions is « 5 X 10 35 erg/s. We used 
the sensitivity for outermost regions because the CXB subtraction 
plays role only in these regions where the surface brightness of the 
source emission is low. This point source sensitivity results in the 
resolved CXB fraction of « 50%, according to Moretti et al. (2003). 
The Galactic and cosmic backgrounds were subtracted from the vi- 
gnetting corrected images and profiles. 

2 http ://cxc . harvard . edu/contrib/maxim/stowed/ 
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Table 1. The list of Chandra observations used for the analysis. 



Obs ID 


^original 




Instrument 


Date 




12 ks 


O.Z 


APT^ T 
A\ v l o - 1 


1 QQQ fVt 1 ^ 

17 77 ULl ID 


305 


4.2 ks 


4.0 ks 


ACIS-I 


1 QQQ TVr 1 1 

i yyy ljcl- i i 


306 


4.2 ks 


4.1 ks 


ACIS-I 


1 QQQ Dp r 97 


307 


4.2 ks 


3.1 ks 


ACIS-I 


9000 Tan 9Q 
Z/UUU J all Z,7 


308 


4.1 ks 


3.7 ks 


ACIS-I 


9000 Ff*h 1 6 


311 


5.0 ks 


3.9 ks 


ACIS-I 


9000 Tnl 9Q 

Z^VjVjVj J 111 Z,7 


312 


4.7 ks 


3.8 ks 


ACIS-I 


9000 Ana 97 


1575 


38.2 ks 


38.2 ks 


ACIS-S 


9001 Ort 05 


1577 


5.0 ks 


4.9 ks 


ACIS-I 


2001 Aug 31 


1583 


5.0 ks 


4.1 ks 


ACIS-I 


2001 Jun 10 


1585 


5.0 ks 


4.1 ks 


ACIS-I 


2001 Nov 19 


2895 


4.9 ks 


3.2 ks 


ACIS-I 


2001 Dec 07 


2896 


5.0 ks 


3.7 ks 


ACIS-I 


2002 Feb 06 


2897 


5.0 ks 


4.1 ks 


ACIS-I 


2002 Jan 08 


2898 


5.0 ks 


3.2 ks 


ACIS-I 


2002 Jun 02 


4360 


5.0 ks 


3.4 ks 


ACIS-I 


2002 Aug 11 


4678 


4.9 ks 


2.7 ks 


ACIS-I 


2003 Nov 09 


4679 


4.8 ks 


2.7 ks 


ACIS-I 


2003 Nov 26 


4680 


5.2 ks 


3.2 ks 


ACIS-I 


2003 Dec 27 


4681 


5.1 ks 


3.3 ks 


ACIS-I 


2004 Jan 31 


4682 


4.9 ks 


1.2 ks 


ACIS-I 


2004 May 23 



In order to study the diffuse emission we need to exclude 
contribution of bright LMXBs. According to the luminosity func- 
tions of LMXBs and faint compact sources associated with old 
population (Gilfanov 2004; Sazonov et al. 2006), the contribu- 
tion of the former is defined by the sources more luminous than 
log L x ~ 35.5 - 36.0. Below this threshold active binaries and cata- 
clysmic variables are the dominating X-ray sources. The individual 
Chandra observations are too short (typically, ~ 4 ksec, see Table 1) 
to achieve this sensitivity, therefore we ran point source detection 
on the combined image with total exposure of T eS = 112.6 ks. To 
combine the data, each observation was projected into the coordi- 
nate system of OBS-ID 303 and the attitude corrections from Voss 
& Gilfanov (2007) were applied in order to better co- align individ- 
ual event lists. To detect point sources we ran CIAO wavdetect 
tool in the 0.5-8 keV band. Some parameters were changed from 
the default values to fit our aims. The scales on which we were 
looking for sources were the V2 -series from 1.0 to 8.0. To mini- 
mize the contribution of residual counts from point sources to dif- 
fuse emission we increased the value of the sigma parameter to 4; 
this parameters describes the size of elliptical source detection re- 
gions in standard deviations assuming a 2D Gaussian distribution 
of source counts. With these adjustment we obtained larger source 
regions than usual, including larger fraction of source counts. To 
increase sensitivity we did not filter out flare containing time inter- 
vals for point source detection. The resulting source list consisted 
totally of 238 sources in the investigated area; the detected sources 
are in good agreement with results of Voss & Gilfanov (2007). The 
sensitivity limit in the central region was 10 35 erg/s, while in the 
outermost region it deteriorated to 5 X 10 35 erg/s. Extracting the 
point spread function with mkpsf for each source we calculated the 
fraction of counts inside the source cell. For most of the sources this 
fraction exceeded 98%, if it was smaller then the source cells were 
enlarged accordingly. The output source cells were used to mask 
out the point sources for further image and spectral analysis. 

We produced exposure maps using a two component spectral 
model consisting of an optically-thin thermal plasma emission with 
temperature of 0.30 keV and a powerlaw component with a slope of 



Table 2. The list of XMM-Newton observations used for the analysis. 



Obs ID 


^original 




Date 


0109270101 


62.5 ks 


16.0 ks 


2001 Jun 29 


0112570101 


61.1 ks 


52.6 ks 


2002 Jan 6 


0112570401 


31.0 ks 


25.0 ks 


2000 Jun 25 


0202230201 


18.3 ks 


17.8 ks 


2004 Jul 16 


0202230401 


14.6 ks 


9.0 ks 


2004 Jul 18 


0202230501 


21.8 ks 


2.0 ks 


2004 Jul 19 



T = 1.80. The ratio of the normalizations of these two components 
was 3/1. This is the best fit two component model to the spectrum 
of the central 200" region. 

2.2 XMM-Newton 

We analyzed 6 archival XMM-Newton observations taken be- 
tween 25.06.2000 and 19.07.2004, as listed in Table 2. We used the 
data of the European Photon Imaging Camera (EPIC) instruments 
(Struder et al. 2001; Turner et al. 2001). For data reduction we used 
Science Analysis System (SAS) version 7.1. 

In order to exclude the flare contaminated time intervals we 
double filtered the lightcurves using hard-band (E > 10 keV) and 
soft-band (E = 1-5 keV) energy ranges according to Nevalainen et 
al. (2005), using 20% threshold for deviation from the mean count 
rate. The remaining useful exposure time is about 58% of the orig- 
inal value. The data was cleaned from the out of time events using 
the Oot event lists. 

The observations were re-projected into the coordinate sys- 
tem of OBS-ID 01 12570101 and merged together. For point source 
removal we combined the source list obtained from Chandra and 
XMM-Newton observations. In regions which lied outside the field 
of view of Chandra we ran the SAS source detection task to com- 
plement the Chandra source list. The source regions were enlarged 
to account for larger size of the point spread function of XMM- 
Newton mirrors, their size was adjusted to approximately match 
~ 90-98% PSF encircled energy radius depending on the bright- 
ness of the point source. It was not possible to reliably exclude 
point source contribution in the crowded central ~ 100" region, 
therefore it was not used in the following analysis. Exposure maps 
were calculated using the eexpmap command of SAS. In trans- 
forming the counts to flux units we assumed the same spectrum 
as for Chandra. 

The particle background on EPIC CCDs consists of two 
components. The "internal" component is generated in inter- 
actions of cosmic rays with the detector material and is ap- 
proximately uniform across the detector. The second compo- 
nent is due to low energy solar protons, concentrated by the 
mirror systems of the telescopes; it is vignetted by the mir- 
rors response but the vignetting is flatter than that for pho- 
tons. The level of both background components is variable 
(see http://www.star.le.ac.uk/~amr30/BG/BGTable.html for de- 
tails). According to this we performed the background subtraction 
in two steps. At the first step the corners of the CCDs which lie 
outside of the field of view, were used to determine the level of the 
flat internal background. The obtained background level was sub- 
tracted from each observations individually. The combined contri- 
bution of the solar protons component and cosmic background was 
approximately determined from the observations of nearby fields 
without extended sources and subtracted from the final vignetting 
corrected image. This method is not perfect, due to the difference 
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Figure 2. Adaptively smoothed Chandra image in the 0.5 - 1.2 keV band 
overplotted with K-band contours. The point sources were removed and 
their locations were filled with the local background value. Center of the 
M31 is marked with a cross. North is up and east is left. 



in the shape of the vignetting function between solar protons and 
photons and due to variability of the solar protons level. However, 
good agreement between emission spectra obtained from Chandra 
and XMM-Newton data confirms adequate accuracy of this proce- 
dure. 



3 RESULTS 
3.1 Images 

The RGB image in Fig.l presents the XMM-Newton data 
(red) overlaid on the 160 \im Spitzer image (Rieke at al. 2004) 
(blue) and K-band image from 2MASS Large Galaxy Atlas (Jarrett 
et al. 2003) (green). Although the main purpose of this image is 
to show the X-ray data coverage, it crudely illustrates the presence 
of a large population of compact sources as well as of the extended 
emission. It also demonstrates the effect of the spiral arms on X-ray 
surface brightness distribution. 

The brightness distribution of the extended emission, after re- 
moval of the point sources, is shown in Fig. 2 along with the con- 
tours of the K-band brightness. The X-ray image was constructed 
from Chandra data in the 0.5-1.2 keV band. The point sources were 
excluded and their locations were filled up with the average lo- 
cal background around the sources. The X-ray image is adaptively 
smoothed. In order to compare the X-ray surface brightness with 
the distribution of the stellar mass we also show K-band contours. 
The surface brightness of the extended emission approximately fol- 
lows the K-band distribution but the image suggests that the agree- 
ment is not perfect and some distortions in the east- west direction 
may be present. In order to investigate this in details we consider 
profiles along the major and minor axes of the galaxy. 




-1000 -500 500 1000 

Major axis offset (arcsec) 




Major axis offset (arcsec) 

Figure 3. Surface brightness distribution along the major axis, in the 
0.5 - 1.2 keV (upper panel) and 2-7 keV (lower panel) energy bands, back- 
ground subtracted. The filled symbols (blue in the color version of the plot) 
show the Chandra data, open symbols (red) show the XMM-Newton data 
and the solid histogram (green) without symbols is the normalized K-band 
brightness. The normalization factors are 4 x 10 27 and 3 x 10 27 erg s _1 
for soft and hard band respectively. The x-coordinate increases from south- 
west to north-east. 



3.2 Surface brightness distribution along the major and 
minor axes 

We studied brightness distribution in two energy bands, 0.5 - 
1 .2 keV and 2-7 keV. Our choice had been motivated by the shape 
of the spectrum of extended emission discussed in the next sec- 
tion. The profiles were constructed along the major and minor axis 
with position angles of 45° and 135° respectively. For each pro- 
file surface brightness was averaged over 500" in the transverse 
direction, corrected for vignetting and the estimated background 
level was subtracted. For XMM-Newton the background level was 
adjusted to achieve better agreement with Chandra profiles. The 
adjusted values were well within the range of the background lev- 
els observed in individual blank- sky observations and differed from 
the average blank- sky level by < 20%. The values of background 
were (4.3, 5.6, 5.2, 6.2) x 10~ 19 erg s -1 arcsec -2 cm -2 for Chan- 
dra (soft, hard) and XMM-Newton (soft, hard) respectively. For all 
distributions we found good agreement between the Chandra and 
XMM-Newton data. We compare X-ray distribution with profiles 
of the K-band emission. As it is well-known, the latter is a good 
stellar mass tracer. The K-band profile was obtained for the same 
regions as the X-ray profiles, in particular the same source regions 
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Figure 4. Same as Fig. 3 but along the minor axis. The normalizations of 
the K-band profile are same as in Fig. 3. The solid histogram with multiple 
peaks (black) shows the distribution of the 160 /urn emission as obtained by 
Spitzer. The x-coordinate increases from south-east to north-west. 



were excluded in computing all profiles. The normalizations for the 
K-band profile are 4x 10 27 and 3 x 10 27 erg s _1 Lq 1 for soft and hard 
band respectively. 

On Fig. 3 we show the distribution along the major axis. In 
the 0.5 -1.2 keV band the profile shows an excess emission in 
the central part of the bulge. At bigger central distances the X-ray 
and K-band light follow each other. In the hard band the X-ray 
surface brightness follows the near-infrared light distribution rather 
well at all central distances, with exception of the shoulder at the 
offset of -300". The excess luminosity of the shoulder above the 
level suggested by the K-band profile is ~ 2 x 10 37 erg/s. Its origin 
is unclear. There is no any easily identifiable feature in the image 
with which it could be associated. It can not be due to residual 
contamination from point sources. Indeed, the excess count rate 
associated with the shoulder is ~ 20% of the total count rate of 
all point sources detected in this region. This is much larger than 
the expected residual contamination from point sources in Chandra 
images, < 2-3%. Good agreement between Chandra and XMM- 
Newton data also excludes the possibility that it is caused by an 
irregularity in the instrumental background. 

The soft band profile along the minor axis (Fig.4) at all off- 
sets exceeds the level suggested by the K-band profile normalized 
according to the X/K ratio from the major axis distribution. More- 
over, the X/K ratio increases significantly outside ~ 500" from the 
center - the X-ray distribution appears to have "wings" extending 
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Figure 5. The X-ray hardness ratio profile along the minor axis, computed 
using the energy bands from Fig.4. The Chandra and XMM-Newton data 
are presented by the same symbols as in Fig.4. The solid histogram (black) 
shows the distribution of the 160 jim emission obtained by Spitzer. The 
x-coordinate increases from south-east to north-west. 

out to ~ 900" or more. Note that exact shape of the surface bright- 
ness distribution at large offsets from the center depends on the 
adopted blank-sky level. The latter can not be directly determined 
from the currently available data, due to its limited field of view. 
For this reason we used the average CXB level and corrected it for 
the fraction of resolved background sources, as described in sec- 
tion 2.1. The obtained value, ~ (4 - 5) x 10~ 19 erg/s/cm 2 , is com- 
parable to the remaining (background subtracted) flux as can be 
seen in Fig.4. Therefore the extend of the X-ray emission at large 
off-center angles can not be unambiguously constrained from the 
present data. In order to eliminate this uncertainty, more extensive 
Chandra observations, including large offset angles are needed. We 
note however, that the existence of unresolved emission at large off- 
sets is independently confirmed by the east- west asymmetry of the 
shadow cast by the 10-kpc star-forming ring, as discussed below. 

There is a clear asymmetry between eastern and western 
halves of the profile, the latter being notably suppressed at offset 
values of +300" and +600". The origin of this asymmetry becomes 
clear after comparison with the 160 iim Spitzer profile. The far- 
infrared profile, plotted in Fig.4 shows several prominent peaks, 
corresponding to the 10-kpc star-forming ring and inner spiral arms 
(Gordon et al. 2006). The surface brightness suppressions in the 
western side of the profile (positive offsets) correspond to the spi- 
ral arm and the 10-kpc ring. On the other hand no significant effect 
of the 10-kpc ring can be seen on the eastern side of the galaxy 
(negative offsets). 

The spatial orientation of M3 1 plays an important role in un- 
derstanding the link between the X-ray and far-infrared distribu- 
tions. We see the galaxy with approximately 77° inclination (Hen- 
derson 1979) and the western side of the galactic disk is closer to 
us (Simien et al. 1978). The effect of this is that we see the western 
side of the bulge through the spiral arms, so the neutral gas and dust 
in the star forming regions cast a shadow on the extended emission 
in the soft band. The eastern side of the disk, on the contrary, is 
located behind the bulge and does not obscure its emission. We es- 
timated the column density of the obscuring material in the spiral 
arms using the observed brightness difference between the eastern 
and western sides and obtained N H ~ 1 - 3 x 10 21 cm -2 . These num- 
bers are compatible with values derived from CO maps (Nieten et 
al. 2006). 
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Figure 6. X-ray spectra of different regions in M31 and of M32: stars (red) 
- inner bulge - the spectrum of central 200" region; filled boxes (green) - 
outer bulge - combined spectrum of two elliptical regions at ~ 700" from 
the center along the major axis; open circles (blue) - spiral arms - two 
elliptical regions at the 10-kpc ring; filled triangles (black) - spectrum of 
M32. The M31 regions are shown in Fig. 1 The inner bulge and M32 spectra 
are obtained by Chandra, outer bulge and spiral arms - by XMM-Newton. 
All spectra are normalized to the K-band flux. 



In the hard band profile along the minor axis we see corre- 
lation between the X-ray and far-infrared emission - the X-ray 
brightness appears to increase at the positions of spiral arms. This 
suggests that spiral arms are sources of harder X-ray emission. To 
further illustrate the impact of spiral arms on the observed X-ray 
brightness we plot in Fig. 5 the hardness ratio along the minor axis 
together with the 160 jdm distribution. This plot confirms the pres- 
ence and significance of the soft emission in the center. The hard- 
ness ratio has clear peaks at the positions of spiral arms, which are 
caused by two effects - obscuration by neutral and molecular gas 
and dust in the soft band and enhanced hard emission associated 
with the spiral arms. 

Based on the X-ray brightness distributions we conclude that 
there are at least three different components in the unresolved X-ray 
emission from the central region of M31: 

(i) broad band component, following the distribution of K-band 
light (i.e. of stellar mass), 

(ii) soft emission, localized in central ~ 500" along the major 
axis of the galaxy and extending out to ~ 900" or more along the 
minor axis, 

(iii) harder emission from the spiral arms and 10-kpc star- 
forming ring. 



3.3 Spectra 

We used the ACIS "blank-sky" files in order to subtract the 
background from Chandra spectra. As before, we renormalized the 
background spectra using the 9.5 - 12 keV band count rates. For 
XMM-Newton we used a combined spectrum of observations of 
nearby fields as a background spectrum and renormalized it using 
the 10 - 12 keV count rate. We found good agreement between 
Chandra and XMM-Newton spectra in all investigated regions. 

In Fig. 6 we show spectra of different regions in M31. The 
spectral extraction regions are depicted in Fig.l. The spectra of 
outer bulge and spiral arms are based on XMM-Newton data, other 



spectra are from Chandra. Also shown is the spectrum of M32. It 
was obtained using Chandra observations OB SID 2017, 2494 and 
5690 with exposure time of r eff = 168.5 ks. The spectrum was ex- 
tracted from elliptical region with the position angle of 40° and with 
90" major and 71" minor axes; the data analysis procedures were 
identical to M31. All spectra are normalized to the same level of 
K-band brightness. The spectral fitting was done in the 0.5 - 7 keV 
band. The element abundances were fixed at Solar value, the hy- 
drogen column density was fixed at the Galactic value. The results 
of spectral fits are summarized in Table 3. 

All spectra shown in Fig. 6 have a prominent power law com- 
ponent with the photon index of approximately r ~ 2 and a soft 
component of varying strength. It is strongest in the inner bulge 
where it by far dominates X-ray emission below 1.2 keV. The best- 
fit temperature of the soft component in the M3 1 spectra obtained 
in the two-component model is in the ~ 0.3 - 0.4 keV range (Table 
3). As can be seen from the table, the simple two-component model 
consisting of a power law and optically-thin thermal plasma emis- 
sion spectrum (MEKAL model in XSPEC) does not adequately de- 
scribe the spectra, especially the spectrum of the inner bulge, hav- 
ing the largest number of counts. The deviations are due to the soft 
band, pointing at the complex shape of the soft component. With 
the second MEKAL component the fit improves significantly, the 
best fit parameters for the inner bulge spectrum are kT\ « 0.2 keV 
and kT 2 * 0.5 keV with^ 2 = 227 for 184 d.o.f. Making the el- 
ement abundance a free parameter improves the fit quality of the 
three-component model further to x 1 - 205 for 183 d.o.f. with 
the best fit abundance of 0.17 of the solar value. We considered 
several other spectral models with free element abundance, they 
all produced sub-solar values in the range of 0.1 - 0.2. We also 
tried to vary element abundances individually and found that the 
model is most sensitive to Ne, Fe and Ni abundances with the best 
fit achieved when the abundances of Ni is a free parameter. The two 
component (vmekal + power law) model requires the Ni abundance 
^3-4 times solar value (abundances of other elements were fixed 
at solar). On the contrary, the models with the free Fe abundance 
give a subsolar best fit values for the latter, ~ 0.6. The fit quality 
improves further with the non-equilibrium thermal emission model 
(vnei model in XSPEC) with the best fit value of n e t ~ 5 • 10 11 
sec/cm 3 and similar dependence on the element abundances as for 
vmekal model. All these modifications do not change the best fit 
temperature significantly. However, none of the models achieves 
acceptable values of x 1 - It is unclear, how much weight should be 
given to these results, as they can be an artifact of the inadequate 
spectral model and insufficient energy resolution of the ACIS -I de- 
tector. Indeed the emission from the inner bulge has a complex 
spectrum composed of several constituents of different temperature 
and ionization state, some of which may be out of the collisional 
ionization equilibrium. 

The outer bulge spectrum has a less prominent soft compo- 
nent, approximately by a factor of 3 weaker than the inner bulge, 
but its temperature, kT « 0.3 keV, is compatible with the in- 
ner bulge value. The spectrum of M32 is similar, although only 
a very faint soft component is present here. Its temperature, kT = 
0.54 + 0.15 keV may be somewhat higher than in M31. All three 
spectra (inner and outer bulge in M31 and M32) nearly perfectly 
match each other above ~ 1.2 keV, after normalization to the K- 
band flux. This is a strong argument in favor of their similar origin. 

The emission from spiral arms clearly stands out. It does not 
have any significant soft component and, most importantly, its nor- 
malization (per unit K-band flux) is by a factor of 4 - 10 higher in 
the hard band than for the other spectra. This difference is smaller in 
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Table 3. Results of the spectral fits in different regions of M31 using Chandra and XMM-Newton observations. PL denotes the power law and MKL is the 
optically-thin thermal plasma emission model. 



Inner region Outer region Spiral arms M32 

Model T/kT ^ 2 /(d.o.f.) r/kT ^ 2 /(d.o.f.) r/kT ^ 2 /(d.o.f.) r/kT ^ 2 /(d.o.f.) 

PL — — — — 2.00 ±0.15 292/192 2.07 ±0.15 113/72 

PL 1.79 ±0.10 2 - 48 !o44 1.65 ±0.21 

MKL 0.36 ±0.01 312/186 0.31 ±0.08 127/101 — — 0.54 ±0.15 83/70 



the soft band, however it is still significantly higher than the spec- 
trum of the outer region. We describe this spectrum with a power- 
law model and we obtain T = 2.0. As with other spectra, the large 
X 1 value indicates more complex spectral shape. 



3.4 Morphology of the soft excess emission 

To explore the spatial distribution of the soft excess emis- 
sion we use two approaches. Firstly we consider the ratio of the 
0.5 -1.2 keV X-ray image to the near-infrared image. The advan- 
tage of the ratio image is that it levels out the component in X-ray 
emission which is proportional to the K-band light. The disadvan- 
tage is that the X/K ratio may become very large at large central 
distances, where the X-ray emission is unrelated to old stellar popu- 
lation and deviates significantly from the K-band light distribution. 
Before producing the ratio image we excluded point sources and 
filled their locations with local background level using the dmf ilth 
tool of CIAO. The X-ray images were adaptively smoothed and the 
near-infrared image was convolved with a gaussian with the width 
comparable to typical smoothing width near the center of X-ray im- 
age. The X-ray images were exposure corrected, the particle back- 
ground and CXB was subtracted. 

Fig.7 shows the ratio images obtained from Chandra and 
XMM-Newton data. They are consistent with each other and illus- 
trate very well the morphology of the soft excess emission. Com- 
parison with the K-band contours shown in the left panel demon- 
strates clearly that the the soft excess emission has nothing to do 
with the old stellar population traced by the near-infrared light. Un- 
like stellar light, it is strongly elongated in the approximate direc- 
tion of the minor axis of the galaxy extending beyond the bound- 
aries of the bulge and projecting onto the 10-kpc starforming ring. 
The images also reveal strong east-west asymmetry, already seen 
in the minor axis distribution (Fig.4). As discussed in the section 
3.2, this asymmetry is caused by absorption of soft X-rays by the 
neutral and molecular gas in the spiral arms on the north-western 
side of the galactic disk. This conclusion is further supported by the 
anti-correlation between 160 iim flux and the soft X-ray brightness, 
with obvious shadows seen at the positions of spiral arms (Fig.7, 
right panel). At the same time, the gas and dust do not have any 
significant effect on the soft X-ray brightness in the south-eastern 
side of the galaxy. 

In the second approach we consider the difference between the 
X-ray and the K-band images, the latter renormalized according to 
the X/K ratio of the old stellar population component. Unlike the 
ratio image it shows the true brightness distribution of the soft ex- 
cess emission. For this purpose we used the same 0.5 - 1.2 keV 
Chandra image as for the ratio image; the normalization of K-band 
image was 4 x 10 27 erg s _1 L" 1 . Fig.8 shows the result. The im- 



age generally confirms the overall morphology of the soft emission 
revealed by the ratio image. 

The images show significant sub- structures in the inner ~ 
100" - 150" with the angular scale of about 30". The origin of 
this substructures is not clear, although some correlation with the 
position of the peaks on the 160 fim image suggests that at least 
some of them may be caused by absorption. These sub- structures 
may deserve a special study which is beyond the scope of this pa- 
per. 



3.5 L X IL K ratios 

We calculate ratios of X-ray to K-band luminosity for the same 
spatial regions as used for the spectral analysis, characterizing the 
inner and outer bulge and spiral arms. As suggested by the spec- 
tra the natural boundary between the soft and hard energy bands is 
1.2 keV. However, in order to facilitate comparison with previous 
studies we computed the ratios in the 0.5-2 keV and 2- 10 keV en- 
ergy bands. The X-ray luminosities for these bands were computed 
using the best fit spectral models. The errors for the luminosity and, 
correspondingly, for X/K ratios account for the model normaliza- 
tion error only and do not include uncertainties in the spectral pa- 
rameters or any other systematic uncertainties. In computing the 
hard band luminosity we extrapolated the best fit model outside the 
energy range used for spectral fits, 0.5-7 keV. The K-band lumi- 
nosities were calculated in each region using the 2MASS image. 

The L X /L K ratios are presented in Table 4. These numbers can 
be transformed to X-ray-to-mass ratios dividing them with the K- 
band mass-to-light ratio, which is of the order of unity. Its value can 
be computed using the B - V « 0.95 colour index from Walterbos 
& Kennicutt (1987) and applying the relation between M+/L K and 
the B-V colour (Bell & de Jong 2001), which gives MJL K * 0.85. 
A close values has been derived by Kent (1992) based on the dy- 
namical mass measurement, M*/L K ^ 1.1. We compare the X/K 
ratios for M31 with M32 and Milky Way. The M32 numbers were 
obtained from the spectra of section 3.3 in the same way as for 
M31. For the Milky Way we used results from Revnivtsev et al. 
(2006) and Sazonov et al. (2006), and transformed them to L X /L K 
for two values of K-band mass-to-light ratio, 0.7 (Dwek et al. 1995) 
and 1 (Kent 1992). The X/K ratio in Sazonov et al. (2006) is given 
for the 0. 1 - 2.4 keV band, following them we converted it to the 
0.5-2 keV band multiplying by 0.7 which is typical ratio for coro- 
nally active stars (Fleming et al. 1995). The contribution of young 
stars is excluded. Revnivtsev et al. (2006) used data from RXTE in 
the 3-20 keV and calculated the 2 - 10 keV luminosities assuming 
a powerlaw spectrum with r = 2.1. 

In the 2 - 10 keV band the X-ray to K-band ratios for bulge 
regions in M31 and for M32 are compatible with each other and 
appear to exceed slightly the Milky Way value. This may be ex- 
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Figure 7. Ratio of 0.5 -1.2 keV image to the the K-band image obtained from Chandra (left) and XMM-Newton (right) data. After the removal of the point 
sources the X-ray images were adaptively smoothed, the K-band image was smoothed correspondingly. The contours in the left and right panels show K-band 
and 160 /urn brightness distributions respectively. The center of M31 is marked with a cross. North is up and east is left. 




Figure 8. The spatial distribution of the soft component. The difference 
between Chandra image in the 0.5 - 1.2 keV band and the K-band image 
normalized according to the X/K ratio of the outer bulge. The contours show 
the intensity levels of 160 /mi Spitzer image. The center of M31 is marked 
with a cross. North is up and east is left. 



M3 1 . Also, in interpreting the soft band ratios the interstellar ab- 
sorption should be taken into account. Indeed, Sazonov et al. (2006) 
studied sources in the Solar neighborhood where the absorption is 
insignificant. The M31 and M32 data, on the contrary, is subject to 
the galactic absorption, with values of 6.1 x 10 20 and 6.3 x 10 20 cm -2 
respectively (Dickey & Lockman 1990). In addition, M31 has spa- 
tially variable internal absorption reaching fewxlO 21 cm -2 in the 
spiral arm regions (Nieten et al. 2006). The total NH for the outer 
bulge region in M31 is probably in the range (7 - 10) x 10 20 cm -2 , 
giving the absorption corrected value of the soft band X/K ratio 
L X /L K * (8 - 10) x 10 27 erg s" 1 L^j. For M32 we obtain after ab- 
sorption correction L X /L K « 4.3 x 10 27 erg s _1 Lq 1 . Taking into 
account the uncertainties quoted by Sazonov et al. (2006) and un- 
certainties in the K-band mass-to-light ratio the X/K ratio for Solar 
neighborhood is formally compatible with the absorption corrected 
value for M31 and probably slightly higher than for M32. It is un- 
clear, how much weight should be given to this discrepancy due to 
the number and amplitude of uncertainties involved. The absorp- 
tion corrected value for the inner bulge of M3 1 is definitely larger, 
L X /L K » (12 - 15) x 10 27 erg s" 1 L^. 

The X/K ratios in the spiral arm region are significantly higher. 
The additional internal absorption can be upto few xlO 21 cm -2 , re- 
sulting in the absorption correction factor in the soft band as large 
as « 4. The final X/K ratios in both bands are approximately an 
order of magnitude larger than for the Milky Way. 



plained by the difference in the spectral shape assumed in the 
flux calculation. Assuming the same power law with T = 2.1 as 
used by Revnivtsev et al. (2006), we obtained for M31 L X /L K = 
(3.0 + 0.1) x 10 27 erg s" 1 L^f, in agreement with the Milky Way 
value. The good agreement of the X-ray to K-band ratios suggests 
similar origin of the 2-10 keV emission in all three galaxies, as 
discussed in more detail in the following section. 

In the soft band, the X/K ratios for M31 are systematically 
larger than M32 ones. This is to be expected from the spectra 
(Fig. 6) and is due to the presence of the soft excess emission in 



4 DISCUSSION 

4.1 Faint compact sources 

We showed that a broad-band emission component exists in 
M31 which (i) follows the K-band light distribution and (ii) its 
2-10 keV X/K ratio is identical to that of M32 and Milky Way. 
These suggest beyond reasonable doubt that this component has 
similar origin to the Galactic Ridge X-ray emission of the Milky 
Way (Revnivtsev et al. 2006). Namely, it is associated with the old 
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Table 4. X-ray to K-band luminosity ratios in different regions in M31 
and in M32 compared to the Milky Way. All values are given in units of 
10 27 , 



erg s" 1 L" 1 . 





0.5 - 2 keV a 


2- lOkeV 


M31 inner bulge 


9.4 ±0.1 


4.5 ±0.1 


M31 outer bulge 


6.2 ±0.1 


3.6 ±0.2 


M31 spiral arms 


22.1 ±0.3 


31.4 ±0.6 


M32 total 


3.5 ±0.1 


4.1 ±0.2 


Galaxy (Sazonov et al. 2006) b 


7.3 - 10.5 (±3) 


2.2-3.1 (±0.8) 


Galaxy (Revnivtsev et al. 2006) b 




2.2-3.1 (±0.5) 



a The X/K ratios for M31 and M32 were not corrected for absorption. See 
the text for the absorption corrected values. 

b The X/M* ratios were converted to X/K ratios assuming two values of the 
K-band mass-to-light ratio, 0.7 (first number) and 1.0; the uncertainty given 
in the parenthesis corresponds to M*/Lk = 1 



stellar population and is a superposition of a large number of weak 
sources of stellar type, the main contributors in the 2-10 keV band 
being cataclysmic variables and coronally active binaries (Sazonov 
et al. 2006). 

The X/K ratios are compatible between all three galaxies in 
the 2-10 keV band, but in the soft band they differ in M31 and 
M32 (Table 4). In the inner bulge this difference is clearly due to the 
contribution of the soft emission from the ionized gas, as discussed 
below. Although we chose the outer bulge region as far as possible 
along the major axis of the galaxy, where the soft X-ray brightness 
follows the K-band profile, some residual contribution from the gas 
can not be entirely excluded. Therefore it is not clear if this differ- 
ence reflects genuine difference between properties and/or content 
of X-ray emitting stellar populations in these two galaxies. Due 
to large systematic uncertainties, the Milky Way value is formally 
compatible with both galaxies and any quantitative comparison is 
inconclusive at this point. 



4.2 Ionized gas 

There are several arguments which suggest that the excess soft 
component has a non- stellar origin. The most important is the mor- 
phology of the excess emission, namely the striking difference from 
the distribution of the near-infrared light (Fig.7). As no significant 
color gradients are observed in the bulge of M31 (Walterbos & Ken- 
nicutt 1988), its stellar content must be sufficiently uniform and 
can not give rise to the observed non-uniformities in the X/K ra- 
tio. Enhanced X-ray to K-band ratio could be explained if a notable 
young population was present in the bulge, which is also not the 
case (Stephens et al. 2003). With the stellar origin excluded, it is 
plausible that the soft excess emission is of truly diffuse nature and 
originates from ionized gas of ~ 0.3 - 0.4 keV temperature. 

To study the physical properties of the ionized gas we use a 
rectangular region on the south-eastern side of the galaxy, to avoid 
the complications due to attenuation by the spiral arms on the north- 
western side. The size of the region is 8' along the major axis and 
10' along the minor axis of the galaxy. From the spectral fit, the 
total X-ray luminosity of the soft component in the 0.5 - 2 keV 
energy range is « 10 38 erg/s, after the absorption and bolometric 
correction the total bolometric luminosity is L bo i « 2.3 X 10 38 erg/s 
assuming Galactic column density. From the emission measure of 
the gas, J n e n H dV = 6.3 x 10 60 cm -3 , we estimate that the mass 
of the gas in the studied volume is of the order of ~ 10 6 M©. We 



note, that assuming the symmetry between the north-western and 
south-eastern sides of the galaxy, the total quantities for the mass 
and bolometric luminosity are twice the quoted values. The average 
number density is about n e ~ 7 X 10" 3 cm -3 . The cooling time of 
the gas is t coo i = (3kT)/(n e A(T)) ~ 250 million years. We also ap- 
plied a two temperature model to the soft emission and obtained for 
the two components: kT\ ~ 0.25 keV, kT 2 -0.6 keV, bolometric 
luminosities of ~ 2 x 10 38 and ~ 10 38 erg/s, total masses of ~ 10 6 
and ~ 0.6 x 10 6 M© and cooling times of ~ 200 and ~ 800 million 
years. We mention that the above computed values strongly depend 
on the applied spectral model. 

The morphology of the gas indicates that it is not in the hydro- 
static equilibrium in the gravitational potential of the galaxy. It sug- 
gests rather, that gas is outflowing from the bulge in the direction 
perpendicular to the galactic disk (see also Li & Wang 2007). The 
mass and energy budget of the outflow can be maintained by the 
mass loss from the evolved stars and Type la supernovae. Knapp 
et al. (1992) estimated the mass loss rate from evolved stars for 
elliptical galaxies ~ 0.0021 L k /L k q M© Gyr -1 . This rate can be 
applicable to the bulge of M3 1 as the stellar populations are similar. 
The K-band luminosity of this region is L K = 1.4 X 10 10 L K ©. The 
estimated total mass loss rate is « 0.03 M©/yr. The stellar yields 
produce the total amount of the observed gas on a timescale of ~ 35 
million years which is shorter than the cooling timescale of the gas. 

To estimate the energy input from type la supernovae we use 
results of Mannucci et al. (2005) who give the supernova rate of 
N S n ia = 0.035+°;°} 3 SNu for E and SO galaxies, where 1 SNu = 
1 SN/10 10 L k q per century. Assuming that one supernova releases 
£sn ia = 10 51 ergs into the interstellar medium, we obtain about 
1.5 X 10 40 erg/s energy that goes into ISM. The minimal energy 
required to lift the gas in the gravitational potential of the galaxy 
can be calculated from E Uft = 12M+o* (David et al. 2006). With 
<r+ = 156 + 23 km/s (Lawrie 1983) we obtain ~ 3.3 x 10 39 erg/s. 
If the gas is heated to the observed temperature by supernovae, it 
requires ~ 10 39 erg/s. This estimates indicate that the energy input 
from supernovae is approximately ~ 3 - 4 times larger than the 
minimal energy required to drive a galactic wind from the galaxy, 
similar to the result of David et al. (2006) for low-luminosity ellip- 
ticals. 

Type Ia supernovae will also contribute to the chemical enrich- 
ment of ISM with iron-peak elements. Typically 0.7 M© of iron is 
provided by each SN Ia event (Nomoto et al. 1984; Shigeyama et 
al. 1992; Iwamoto et al. 1999), that gives about 3.4 x 10" 4 M© of 
iron per year. Assuming complete mixing of the supernova ejecta 
with the stellar wind material we would expect the iron abundance 
in the hot ISM « 1.1 x 10~ 2 by mass, which exceeds the Solar value 
of 1.9 x 10" 3 (Anders & Grevesse 1989) by a factor of ~ 6. The ob- 
served spectra are inconsistent with high iron abundance in simple 
one- or two-component thermal models, but these models can not 
adequately describe the spectra anyway, therefore this result can 
not be used as a conclusive argument. On the other hand, the dis- 
crepancy between high predicted and low observed abundance of 
iron is a well-known problem for elliptical galaxies, where type Ia 
supernovae also plays important role in thermal and chemical evo- 
lution of the ISM; it has been addressed in a number of studies, e.g. 
Brighenti & Mathews (2005). We also note, that although the iron 
is the primary element by mass in the type Ia supernovae ejecta, the 
ISM will be also enriched by other elements, most significantly by 
nickel. The detailed analysis of this problem is beyond the scope of 
the present paper. 

Assuming that the gas leaves the galaxy in a steady state wind 
along the axis perpendicular to the plane of the disk, the outflow 
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speed v w can be calculated from M* = nr 2 p gas v w where p gas is the 
average gas density estimated above and r is the radius of the base 
of the imaginary cylinder filled with the outflowing gas. This cal- 
culation gives v w ~ 60 km/s which is smaller by a factor of few 
than the local sound speed. Such a slow sub- sonic motion of the 
gas can not explain the observed elongated shape of the gas distri- 
bution, which may be related to the the magnetic fields and galactic 
rotation. We can use the fact that the shadow from the 10-kpc star- 
forming ring is present to estimate the extent of the gas along the 
axis perpendicular to the galactic disk. The angular distance of the 
shadow from the center is ~ 600" - 700", giving the "vertical" 
extent of the gas of > 2.5 kpc. 



4.3 Spiral arms 

The -10 times higher X-ray to K-band ratios observed in the 
spiral arms (Table 4) and their different emission spectrum (Fig. 6) 
suggest that X-ray emission from spiral arms has different nature 
than the bulge. As spiral arms are associated with star-formation, 
an obvious candidate is X-ray emission from young stellar objects 
(protostars and pre-main sequence stars) and young stars, which 
are well-known sources of X-ray radiation (Koyama et al. 1996; 
Carkner et al. 1998; Giardino et al. 2007). 

As X-ray emission from the spiral arms is associated with 
young objects, it is natural to characterize it with L X /SFR ratio. 
We compute this ratio for the regions used in spectral analysis. The 
FIR flux was determined from the 160 fim Spitzer image, 290 Jy. 
In computing this value we subtracted the blank- sky background of 
nearby fields. To convert it to SFR we used results of IR spectral fits 
from Gordon et al. (2006), which gave SFR = 9.5 x 10" 5 F^/Jy 
M©/yr for the M31 distance. Thus we obtained the star-formation 
rate of 0.028 M©/yr in the region used for the analysis. The X-ray 
luminosity in the same region is 3 X 10 37 erg/s in the 2 - 10 keV 
band. After subtracting the X-ray emission due to the old stel- 
lar population we obtain 2.6 x 10 37 erg/s energy. From this we 
can compute L X /SFR * 9.4 x 10 38 (erg/s)/(M Q /yr). This value 
is ~ 1/3 of the total L X /SFR arising from HMXBs which is 
2.5 x 10 39 (erg/s)/(M /yr) (Grimm et al. 2003; Shtykovskiy & 
Gilfanov 2005). 



5 CONCLUSION 

We investigated the origin of unresolved X-ray emission from 
M31 using Chandra and XMM-Newton data. We demonstrated that 
it consists of three different components: 

(i) Broad-band emission associated with old population, similar 
to the Galactic ridge emission in the Milky Way. It is a a combined 
emission of a large number of weak unresolved sources of stellar 
type, the main contribution being from cataclysmic variables and 
active binaries. The surface brightness distribution of this compo- 
nent approximately follows the distribution of K-band light. The 
absorption corrected X-ray to K-band luminosity ratios are com- 
patible with the Milky Way values. The total luminosity of this 
component inside central 800" x 400" is of the order of ~ 3 x 10 38 
erg/s in the 0.5 - 10 keV band. 

(ii) Soft emission localized in the inner bulge of the galaxy along 
its minor axis. This emission is from ionized gas with the tempera- 
ture of the order of ~ 300 eV, although its spectrum can not be ad- 
equately described by a simple one- or two-temperature model of 



optically-thin emission from a gas in a collisional ionization equi- 
librium. The 0.5 - 2 keV luminosity in the central 8' x 20' area is 
~ 2 X 10 38 erg/s, the absorption corrected bolometric luminosity 
is ~ 5 x 10 38 erg/s. The total mass of the gas is ~ 2 x 10 6 M©, 
its cooling time ~ 250 Myrs. The surface brightness distribution 
is drastically different from the stellar light distribution, it is sig- 
nificantly elongated along the minor axis of the galaxy. The mor- 
phology of the soft emission suggests that gas outflows along the 
direction perpendicular to the galactic plane. The mass and energy 
budget is maintained by the mass loss by the evolved stars and type 
la supernovae. The "vertical" extent of the gas exceeds > 2.5 kpc. 
These results are in good agreement with those recently obtained 
by Li & Wang (2007). 

(iii) Hard emission from spiral arms. In the 0.5-7 keV band 
this emission has approximate power law shape with r ~ 2. It 
is most likely associated with star-forming regions and is due 
to young stellar objects and young stars. The L X /SFR « 9.4 x 
10 38 (erg/s)/(M©/yr) which is about ~ 1/3 of the contribution of 
HMXBs. 
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